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Abstract 



Using 13.4 fb of data collected with the CLEO detector at the Cornell Elec- 
tron Storage Ring, we have observed 300 events for the two-photon production 
of ground-state pseudo-scalar charmonium in the decay n c — > K^K^ir^. We 
have measured the n c mass to be (2980.4 ± 2.3 (stat) ± 0.6 (sys)) MeV and 
its full width as (27.0 ± 5.8 (stat) ± 1.4 (sys)) MeV. We have determined the 
two-photon partial width of the n c meson to be (7.6 ± 0.8 (stat) ± 0.4 (sys) ± 
2.3 (br)) keV, with the last uncertainty associated with the decay branching 
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In this Letter, we report a study of two-photon production of the ground-state pseudo- 
scalar charmonium, i.e., 77 — > r] c . The two space-like photons are radiated by e + and e~ 
beams, each at an energy of approximately 5.3 GeV. The charmonium spectrum is an ideal 
testing ground for quantum chromodynamics (QCD) calculations, and producing C-even 
charmonium states through 77 fusion provides a clean environment for this purpose. 

The two-photon partial width of the r] c meson can be expressed in next-to-leading order 
(NLO) perturbative QCD (PQCD), in terms of the e + e~ partial width of the J/ip meson, 
as! 

I> 4 |^ foil 2 

= -(1 + 1.96a a /7r) x 1 ric ) i . (1) 



Using the world average value || of Tf e , a value of the strong coupling constant a s evaluated 
at the charm mass scale |I[ of (0.28 ± 0.02), and the assumption that the two IS wave 
functions, are the same at the origin, this relationship predicts TUl = (8.2 ± 0.6) keV. 

The total width of the r] c meson can be assumed to be dominated by its two-gluon 
component, i.e., F^t ~ H^. The ratio of the rates for r\ c — > gg and r\ c — > 77 is an especially 
clean prediction of PQCD because the dependencies of these rates on the wave functions 
and non-perturbative factors are identical in the numerator and denominator. The ratio 
depends only on the coupling constants and has been calculated in NLO M, 

Ltt „ 9a? (l + 4.8a 5 /7r) 

~ 8a 2 (l-3.4a s /7r)' 1 ' 

Using the value of TUl estimated in NLO gives r^t as (28 ± 6) MeV; using instead the 
world average value of one obtains an estimate of r^ t as (26 ± 6) MeV. A calcu- 
lation with fully relativistic decay amplitudes and a sophisticated QCD potential model ]3[] 
predicts r^t ~ 23 MeV. The current world average of r^ t = 13.2i|J MeV disagrees with 
these theoretical expectations. A precise measurement of the full width and two-photon 
partial width of the rj c is important for the verification of these PQCD calculations and 
approximations. 

In the two-photon process e + e~ — > e + e _ 77 — > e + e~r] c , the photon propagators dictate 
that the cross section naturally peaks at low momentum transfer, so the photons are almost 
real ( "on shell" ). The incident leptons are scattered at very low angles and continue traveling 
down the beam pipe undetected. Such "untagged" events typically have low net transverse 
momentum and low visible energy. The production of the 7] c meson in this untagged two- 
photon process was searched for in the K^K^n^ 1 decay mode. 

The data used in this study correspond to an integrated luminosity of 13.4 fb _1 and were 
collected with two configurations (CLEO II [§ and CLEO II. V 0) of the CLEO detector at 
the Cornell Electron Storage Ring (CESR). Approximately one third of the data were taken 
with the CLEO II configuration. The detector components most useful for this study were 
the concentric tracking devices for charged particles, operating in a 1.5T superconducting 
solenoid. For CLEO II, this tracking system consisted of a 6-layer straw tube chamber, 
a 10-layer precision drift chamber, and a 51-layer main drift chamber. The main drift 
chamber also provided measurements of the specific ionization loss, dE/dx, used for particle 
identification. For CLEO II. V, the straw tube chamber was replaced by a 3-layer, double- 
sided silicon vertex detector and the gas in the main drift chamber was changed from a 
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50:50 mixture of argon-ethane to a 60:40 helium-propane mixture. These changes gave rise 
to significant improvements in the momentum and dE/dx resolutions for charged tracks. 
Photons were detected using the high-resolution electromagnetic calorimeter consisting of 
7800 Csl crystals. The Monte Carlo simulation of the CLEO detector response was based 
upon GEANT [[J. Simulated events were processed in the same manner as the data to 
determine the rj c — ► KgK^n^ detection efficiency and the KgK^^ mass resolution at the 
rj c mass. 

The K s vertex was reconstructed from its decay to 7r + 7r~ and was required to be dis- 
placed from the e + e~ interaction point; the amount of this displacement varied with detector 
configuration but was ~1.5 mm. The K s candidate was also required to be within 4 standard 
deviations (a) of the known K s mass [[| ; here o was determined on an event-by-event basis 
from the momenta measurements. Furthermore, the K s momentum vector was required to 
point back to the interaction point. Of the two remaining charged tracks, the K T and 7r 
candidates, the one with lower momentum was typically uniquely identified using the parti- 
cle's specific ionization (dE/dx). This fixed the identity of the only remaining unidentified 
track, because the presence of the Kg dictated that exactly one of these two be a kaon to 
conserve strangeness in the event. 

The background from processes other than two-photon production was suppressed by 
requiring that the r] c candidate have net transverse momentum less than 0.6 GeV/c and 
that visible energy in the event be less than 6 GeV. Also, because the final state had no 
expected energy deposits in the calorimeter from neutral particles, the total calorimeter 
energy in the event not matched to charged tracks was required to be less than 0.6 GeV. 

For the mass measurement only, we restricted ourselves to those events in which the K T 
and 7r daughters of the r] c candidate traversed all layers of the tracking volume. The K s 
daughter pions were not required to satisfy the same criterion, in that the kinematic fitting 
of the K s decay corrected for any possible momentum mis-measurement of the daughter 
tracks. We did not make any such requirements while measuring T^ t and T^, because 
these quantities are relatively insensitive to precise measurements of the track momenta; 
the distribution of candidate invariant masses for the determination of these two quantities 
is shown in Fig. [I]. 

We fitted the background with a power law function (A ■ W^, with the K s K T ii ± 
invariant mass and A a multiplicative constant) and the signal with a spin-0 relativistic 
Breit-Wigner function (describing the natural line shape) convolved with a double Gaussian 
function (to take into account the detector resolution). The parameters for this double 
Gaussian were obtained from a Monte Carlo sample that had the r] c intrinsic width set 
to zero. We performed a simultaneous, binned, maximum-likelihood fit to the invariant 
mass distributions for the CLEO II and CLEO II. V datasets, constraining the physical 
variables M Vc , r^ t , and in the two datasets to be the same. The constraint on r^s, was 
accomplished by requiring the ratio of the fitted yields to be the same as the ratio of the 
integrated luminosities of the two data sets times the efficiencies as determined from our 
simulations. The invariant mass resolution was approximately 9 MeV in CLEO II and 7 
MeV in CLEO II. V. The bin width for fitting was chosen as approximately the average of 
these two resolutions. 

As noted above, two separate sets of fits were performed, one for the measurements of 
r^t and T^ c and another for the determination of the rj c mass. The full width and yield 
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were obtained from the distributions shown in Fig. [T], with the total observed yield being 
N obs = 300 ± 32. The fit to the width gives T^ t = (27.0 ± 5.8) MeV, with the uncertainty 
being only statistical. The two-photon partial width was determined by first correcting for 
the detector efficiency, e, and then dividing by the number of events expected, Aq, for a 
two-photon partial width of 1 keV: 



H c 7 = N obs /{e ■ N x 



(3) 



The quantity Aq was determined using 



N\ = £ • B Vc • Bk s ■ (Je+e-^e+t 



"Id 



(4) 

e + e~rj c was obtained from 
7§ and setting Ffc = 1 keV; 
this choice of this value has no effect on the extracted result. Also, B Vc = B(r] c —+ K^K^ir^) 



£ is the integrated luminosity. The cross section for e + e 
Monte Carlo simulation, using the formalism of Budnev et al. 
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and Bk s = B(K$ — > with the world average values used. This procedure gives 

(7.6 ± 0.8) keV, with this uncertainty coming from statistics only. 
The mass was obtained from fits to the more restrictive set of events, as described above, 
with a total signal size of 195 ± 24, yielding M Vc = (2980.4 ± 2.3) MeV, the uncertainty 
being statistical only. 



1600500-002 




M(Kg \lCjt ) 

FIG. 1. Results of the simultaneous fit to r/ c candidates in (a) CLEO II and (b) CLEO II. V for 
width and yield measurements, with a combined x 2 /d.o.f. = 226/243. 



Possible sources of systematic uncertainty for the measured mass, full width, and two- 
photon partial width were studied. The results are summarized in Table I, in which the 
individual uncertainties are added in quadrature to obtain the total systematic uncertainty. 

The mass calibration of our detector was checked by measuring the masses of the well 
known Kg, 0, D, and J/ip mesons using decay modes involving only charged tracks. The 
measurements were found to be in good agreement with their respective world averages when 
we limited ourselves to events in which the charged tracks traversed all layers of the tracking 
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volume. The fitted mass including events outside this more restrictive detector volume is 
consistent with the value we quote, but with substantially larger systematic uncertainties. 
Our particle identification procedure did not introduce any significant systematic bias to 
the mass measurement. The measured mass was also insensitive to variations in the signal 
shape used to fit the r\ c resonance. 

The systematic uncertainties in the measurement of the full width were dominated by 
effects due to the mass resolution, the particle identification procedure, and the signal shape 
used in the fit. Our ability to predict the actual mass resolution was tested by studying the 
reconstructed D mesons in Monte Carlo simulation and data; the agreement was found to 
be better than 0.1 MeV. 

The particle identification procedure was unable to distinguish between a charged kaon 
and a charged pion if the track momentum was above 0.8 GeV/c, for which the expected 
ionization losses are nearly equal for the two species. This led to a broadening of the recon- 
structed resonance and was taken into account by the wider Gaussian of the double Gaussian 
resolution function. This effect was limited to less than 5% of the events. We estimated the 
possible uncertainty due to this mis-assignment of particle species by completely removing 
the fraction of events having two possible r\ c candidates and assigning the corresponding 
change in the measured width as the systematic error. 

The accuracy of the fitting method in extracting the Breit-Wigner width of the resonance 
was checked by extracting the r\ c widths from sets of simulation events generated with 
different intrinsic widths. We varied the parameters of the signal shape used to fit the 
f] c resonance within their uncertainties, derived from a comparison of the fit to D meson 
decays in Monte Carlo simulation and data, to estimate the effect on the measured width. 
The measured width from the more restrictive set of events used for mass measurement was 
within 0.2 MeV of the corresponding measurement using the full sample, and no significant 
correlation was found between the measured mass and full width. 

There were several sources of uncertainty for the estimation of the efficiency, which in 
turn affected the measurement of H^, . These were dominated by the uncertainties in the 
tracking and trigger efficiencies. The effect of a possible presence of resonant substructure 
(K*K) in r] c decay was studied and was found to give an insignificant variation in the 
detection efficiency. We estimated the uncertainty in the measured partial width from this 
effect by considering the possibility that all the i] c mesons decay through K*K. Our initial 
investigation showed roughly a third of the i] c — > K^K t tt ± events proceed via If* (1430) 
with no evidence for any fT*(892); a detailed analysis of this substructure is beyond the 
scope of this Letter. 

The selection requirements on total visible energy and unmatched energy clusters were 
shown by simulation to be essentially 100% efficient for our signal process and free of sys- 
tematic bias. Possible bias from the transverse momentum requirement was investigated by 
changing the nature of the form factor in the simulation and shown to also be negligible. 

In our analysis, we investigated the possible effects of interference between the KgK T n ± 
resonant and non-resonant final states. From the distribution of net transverse momentum 
for events in the sidebands of the signal lineshape, we estimated that one third of the 
background events were not of the type 77 — » K^K^n^; these included events with at least 
one missing particle (7r , 7) as well as events of the type e + e~ — > hadrons, e + e~ — > t + t~, 
and 77 — > t + t~. Such events cannot have interference with our signal KgK T ir ± events. 
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Study of the helicity angle distributions in the 7] c rest frame indicated that the sideband 
events predominantly have J = 2 (or higher) while our signal has J = 0. Due to the 
preferential production of states with natural parity in two-photon un-tagged processes, a 
majority of the remaining background events were expected to have the natural spin-parity 
(0 + ) compared to the unnatural spin-parity (CT) of the signal events. The acceptance of 
our detector is symmetric in polar angle and uniform in azimuth, making the interference 
between these states of opposite parity vanish. We did not include any possible effects due 
to interference on the measured mass, full width and two-photon partial width. Further, the 
signal shape showed no distortions and the goodness of fit to the hypothesis that ignored 
interference was very good, as shown in Fig. [L| 

In summary, we have measured the mass, full width, and two-photon partial width of the 
T] c produced in two-photon collisions. The mass measurement of (2980.4 ± 2.3 (stat) ± 0.6 
(sys)) MeV compares well with the world average @ of (2979.8 ± 2.1) MeV. The measured 
total width of (27.0 ± 5.8 (stat) ± 1.4 (sys)) MeV disagrees with the world average || of 
(13.2l|;§) MeV, which consists of measurements with large relative uncertainties (40—100%). 
Our measured width is consistent with theoretical expectations The measured two- 

photon partial width of the r] c meson of (7.6 ± 0.8 (stat) ± 0.4 (sys)) keV agrees well with 
the world average of (7.5l};4) keV and theoretically expected values, and is a significant 
improvement in terms of experimental precision. We use the world average [[J of the 7] c — > 
KKtc branching fraction of (5.5 ± 1.7)%. The uncertainty in due to the uncertainty in 
this branching fraction is ±2.3 keV and is stated separately from the other contributions. 
From the ratio of our measured full width and two-photon partial width, we have extracted 
a s at the charm mass scale to be 0.285 ± 0.025, for which we have added our sources of 
uncertainty in quadrature. We have used the NLO calculation in Eq. 2 to estimate a s , thus 
making the result dependent on renormalization scheme and scale; we have not included 
such theoretical uncertainties in our quoted value. 

Our measurements of T^t and show that PQCD calculations are able to reliably 
predict the ratios of the decay rates of a heavy quarkonium system, where non-perturbative 
effects cancel. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences 
and Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss Na- 
tional Science Foundation, the Texas Advanced Research Program, and the Alexander von 
Humboldt Stiftung. 
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Source of 


M Vc 


tot 


77 


systematic uncertainty 


(MeV) 


(MeV) 


(keV) 


Mass calibration 


0.6 


< 0.1 


< 0.1 


Particle identification 


< 0.1 


1.3 


0.1 


Signal shape 


< 0.1 


0.3 


0.1 


Detector resolution 


< 0.1 


0.3 


< 0.1 


Trigger 






0.2 


Tracking 


- 


- 


0.2 


Resonant substructure 


< 0.1 


< 0.1 


0.2 


Luminosity 






0.1 


Kg selection 


< 0.1 


< 0.1 


0.1 


Event Selection 


< 0.1 


< 0.1 


< 0.1 


Total 


0.6 


1.4 


0.4 



TABLE I. Systematic uncertainties in the three measurements. The overall value is obtained 
by adding the individual contributions in quadrature. 
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